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ABSTRACT 

Due to its centrally bright X-ray morphology and limb brightened radio profile, MSH 11-61A (G290.1-0.8) 
is classified as a mixed morphology supernova remnant (SNR). Hi and CO observations determined that the 
SNR is interacting with molecular clouds found toward the north and southwest regions of the remnant. In 
this paper we report on the detection of 7 -ray emission coincident with MSH 11-6IA, using 70 months of 
data from the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. To investigate the origin 
of this emission, we perform broadband modelling of its non-thermal emission considering both leptonic and 
hadronic cases and concluding that the 7 -ray emission is most likely hadronic in nature. Additionally we 
present our analysis of a 111 ks archival Suzaku observation of this remnant. Our investigation shows that the 
X-ray emission from MSH 11-61A arises from shock-heated ejecta with the bulk of the X-ray emission arising 
from a recombining plasma, while the emission towards the east arises from an ionising plasma. 

Subject headings: ISM: individual (MSH 11-61A, G290.1-0.8) — cosmic rays — X-rays: ISM — gamma 
rays: ISM — supernovae: general — ISM: supernova remnants 


1. INTRODUCTION 

Since the launch of the Large Area Telescope (LAT) on¬ 
board the Fermi Gamma-ray Space Telescope, its improved 
sensitivity and resolution in the MeV-GeV energy range has 
lead to a number of supernova remnants (SNRs) being de¬ 
tected in GeV 7 -rays. The shock-front of an SNR is expected 
to be able to accelerate cosmic rays (CRs) efficiently, pro¬ 
ducing non-thermal X-ray and 7 -ray emission. As 7 -rays can 
arise from leptonic processes such as inverse Compton (IC) 
scattering and non-thermal bremsstrahlung from high-energy 
electrons, or from hadronic emission arising from the decay 
of a neutral pion (produced in a proton-proton interaction) 
into two photons, a means of distinguishing between these 
two mechanisms is crucial for our understanding of the origin 
of this observed emission. Thermal and non-thermal emission 
from SNRs have provided increasing support in favour of CRs 
bein g accelerated at the shock front of the r emnant (e.g. Ty¬ 
cho: JWarren et al.lHool RX J1713.7-39 46: lUchivama et all 
120071 W44. MSH 17-39 and G337.7-0. 1 : fCastro et alJUoTlf 
SNRs known to be interacting with dense molecular clouds 
(MCs) are ideal, indirect laboratories that one can use to de¬ 
tect and analyse 7 -rays arising from accelerated protons. The 
interaction of the SNR’s shockwave with dense molecular ma¬ 
terial is often inferred from the detection of one or more OH 
(1720 MHz) masers, but enhancement of excitation line ra¬ 
tios such as 7 = 2 —> 1/7 = 1 —>• 0, broadenings and asymme¬ 
tries in molecular line features or morphological alignment 
of molecular features with SNR features can also allow one 
to determine SNR/MC interaction (see lSlane et al.112015l and 
references therein). 

Som e of the firs t SNRs detected by the Fermi-L AT (e.g. 
W44: lAbdoet all 1201 Obi lAckermann etal] I 20 TI IC 443: 
(lAckerma nn et alj|2013[): 3C391: ICastro & Slanell 2010 l and 
W49B: lAbdo et al.ll 2010 ah are part of a unique class called 
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Mixed-Morphology (MM) SNRs. Some of these SNRs are 
known to be interacting with MCs. These SNRs are char¬ 
acterised by their centrally peaked X-ray morphology which 
is thermal in n ature, while their radio profiles are limb- 
brightened dRho & Petrel 11998 ). The evolutionary sequence 
leading to these unusual X-ray properties are not well un¬ 
derstood and the morphology and characteristics of these 
SNRs are difficult to explain using standard SNR evolution 
models. There are two main models that are invoked in 
the literature to exp lain their characteristics. One possible 
model (White & Long! 11991 ) assumes that in the vicinity of 
the supernova explosion there are many small, dense, cold 
cloudlets. These cloudlets are small enough that they do not 
affect the passage of the shock, and are sufficiently dense that 
they are neither blown apart nor swept up. Once the shock 
has passed, the cloudlets slowly evaporate, filling the interior 
of the SNR with a relatively dense gas that emits in X-rays. 
Another possible scenario is that thermal conduction results in 
the transport of heat and material to the center of the remnant, 
increasing the central density of the remnant, and smoothing 
the temperature gradient behind the shock (Cox et ~aklfl999l) . 

The ionisation state of a thermal plasma in an SNR can 
be characterised by its ionisation temperature ( Tz ) which de¬ 
scribes the extent to which the ions are stripped of their elec¬ 
trons and its electron temperature (T e ) which describes the ki¬ 
netic energy of the electrons. The thermal plasmas of SNRs 
have been thought to be either underionised, where Tz < T e 
or in collisional ionisation equilibrium Tz = T e . Recent ob¬ 
servations by the Suzaku satellite have confirmed earlier sug¬ 
gestions based on ASCA data (Kawasaki et al.ll2005 ). that the 
thermal plasma in som e MM SNRs is overionised (recom - 
bining) (e.g. 3C391: ISato et all 120141 lErgin et all 1201 4) . 
Recombining plasmas have ionisation temperatures that are 
higher than the electron temperatures and require rapid cool¬ 
ing of electrons either by thermal conduction (Kawasaki et all 
| 2002l), adiabatic exp ansion via rarefraction and recombination 
(Itoh & Masail 19891 or the interaction with dense cavity walls 
or molecular clouds (Dwarkadasll2005l) . 

MSH 11-61A (G290.1-0.8) is a Galactic MM SNR that 
is known to be interacting with MCs. It was first discov- 
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ered by iMills et ahl (119611) using the Sydney 3.5m cross¬ 
type radi o telescope. It was first identified as an SNR by 
iKestevenl (1 19681) and later classified as a shell-type SNR 
with a complex internal structure and ear-like protrusions to¬ 
wards the northwest and southeast using the Molongo Ob¬ 
servatory Syn thesis Telescope (MOST) at multiple different 
wavelengths dKesteven & Caswelll 119871 iMilne et ahl 119891 
IWhiteoak & Green! 1 19961) . It has an angular size of 19’ x 
11 ’ and a radi o-co ntin uum spectral index of a = -0.33 + / — 
0.07 dRevnoso et al.l 12006 ). Radio continuum observations 
using the Australia Telescope Compact Array (ATCA) by 
[Reynoso et al.l ( 2006 ) showed filamentary emission with lit¬ 
tle shell structure, while the northern and southern edges of 
the remnant show evidence that the shock front could be in¬ 
teracting with a plane parallel density gradient. Using NAN- 
TEN CO images of MSH 1 l- 6 /A. [Filinovic et~al J (12005!) de¬ 
termine that the SNR is associated with a MC towards the 
south-west and northern rim of the remnant. Hi observations 
using the Southern Galactic Plane Survey find that the molec¬ 
ular cloud is found at a local standard of rest velocity of ~ 13 
km s -1 (iMcClure-Griffiths et al]|2005l) . 

MS H 11 -6 / A was first detected in X-rays by [Seward 
(119901) using a 10.9 ks Einstein Observatory observation. The 
0.3 - 4.5 keV Imaging Proportional Counter (IPC) image of 
the remnant showed that the X-ray emission is peaked to¬ 
wards the center of the remnant. Using a 40 ks Advanced 
Satellite for C osm olo gy and Astrophysics ( ASCA ) GIS ob¬ 
servation, ISlane et all (12002 ) were able to determine that the 
central X-ray emission is thermal in nature, classifying the 
remnant as a MM SNR. They modelled the X-ray emis¬ 
sion using the cloudy ISM model by [White & Long! (119911) 
and derived an intercloud medium density of ~ 0.05 - 0.40 
cm ~ 3 and an age of 10 - 20 kyr. Using XMM-Newton, 
[Garcia et al.l (120 12l) analysed five regions along the axes of 
the remnant using an absorbed plane parallel non-equilibrium 
ionisation (VPSHOCK) model and found that the physical 
conditions across the remnant are not homogeneous, with 
variatio n in ionisation state , tempe rature and elemental abun¬ 
dances. iKamitsukasa et al.l (1201 5l) analysed Suzaku data and 
found that in the center and in the northwest of the remnant 
the plasma is recombining, while everywhere else it is ionis¬ 
ing. 

The distance to MSH 11-61A has been measured by a num¬ 
ber of different meth ods. Hi measurements taken by the 64- 
m Parkes telescope dDickell 1 1973h iGoss et elL| 1 1 972ft gives a 
lower l imit o f ~_3.5 kpc to the remnant. Ha measurements by 
iRosado et al.l d 1996b using a Fabry-Perot interferometer im¬ 
plied a distance of 6.9 kpc assuming a Vlsr o f +12 kmjf " 1 
for the SNR. Using CO measurements. IRevnoso et a ll (120061) 
derive d a distance of 7 - 8 kpc assum i ng the iBrand & Blita 
(119931) rotation curve. IRevnoso et al.l (120061) derived a dis¬ 
tance of 7 ± lkpc using Hi absorption measurements from 
ATCA, combined with data from the Southern Galactic Plane 
Survey, while lSlane et al.l (120021) estimated a distance of 8 -11 
kpc by modelling the thermal X-ray emission of the remnant 
as detected by ASCA. We use 7 kpc throughout this paper. 

There are three pulsar s close to the position of 
MSH 11-61A. IKaspi et all (fl997h discovered the young 
(spin-down age, r = 63 kyr), energetic pulsar PSR J1105- 
6107 (J 1105) which is located approximately 25’ away from 
the remnant. It has a spin-down luminosity of 2.5 x 10 36 erg 
s -1 and overlaps the position of the EGRET 7 -ray source 3EG 
J1103-6106. It was also detected via periodicity searches in 


GeV 7 -rays by the Fermi-LAT satellite (lAbdo et al .1120131) . It 
has a dispersion measure of 271 cm -3 pc which implies a dis¬ 
tance o f ~7 kpc using the stan da rd Galactic e lectron density 
model (ICordes & Laziol [20021) . IKaspi et ah] (119971) consid¬ 
ered the scenario that this pulsar is associated with the rem¬ 
nant and determined from proper motion measurements that 
it would need to be travelling with a transverse velocity of ~ 
650 km s -1 to have reached its current position, assuming a 
distance of 7 kpc to the pulsar and r = 63 kyr. This is much 
larger than the average pulsar transverse velocity but much 
less than what has_been sugg ested for other pulsar-SNR asso¬ 
ciations (e.g JCaraveolll993i) . leading the authors to conclude 
that association is possibl e. Usi ng the ASCA X-ray charac¬ 
teristics of MSH 11-61A. [Slane et ahl (120021) concluded that 
MSH 11-61A and J1105 are not associated under the as¬ 
sumption that the SNR evolved via thermal conduction or a 
cloudy ISM. These two models imply a transverse velocity 
of ~ 4.5 x 10 3 km s -1 and ~ 5.3 x 10 3 km s -1 respectively, 
which is much larger than the mean velocity (~ 310 km s -1 ) 
of young pulsars oHobbs et al.ll2005l). The two other pulsars, 
PSR Jl 103-6025 dKramer et al.l 120031) and PSR J1104-6103 
(IKaspi et al .11 19961) are not associated with the remnant as they 
have characteristic ages much larger than 1 Myr, which is far 
greater than the expected lifetime of an SNR. Also nearby is 
the INTEGRAL source ICG Jl 1014-6103, which is a neutron 
star trav ellin g at a velocity exceeding 1000 km/s, which the 
iPavan et aD (120141) associate with MSH 11-61A. 

Using 70 months of Fermi-LAT data, we analyse the GeV 
7 -ray emission coincident with MSH 11-61A and investigate 
the nature of this emission using broadband modelling. In 
addition, we analyse archival Suzaku data and report on the 
spatial and spectral properties of the X-ray emission of this 
remnant. In Section 2 we describe how the Fermi-LAT data 
are analysed and present the results of this analysis. In Section 
3 we present our spatial and spectral analysis of the Suzaku 
observation of MSH 11-61A, while in Section 4 and 5 we 
discuss the implications of our results. 

2. OBSERVATIONS AND ANALYSIS OF MSH 11-61A 

We analysed ~70 months of reprocessed data collected 
by the Fermi-LAT from 4th August 2008 to 16th June 
2014. We selected data within a radius of 20° centered on 
MSH 11-6 JA. We used the “P7REP_SOURCE_V15” instru¬ 
ment response function (IRF) which is based on the same in¬ 
flight event analysis and selection criteria that was used to 
generate the previous “P ASS7_V6” IRFs (details described in 
lAckermann et al.l[2012 ). Due to the improved reconstruction 
of the calorimeter position as well as a 1 % per year correc¬ 
tion for the degradation of the light yield of the calorimeter, 
the new IRFs significantly improve the point spread func¬ 
tion (PSF) of the LAT for energy > 5 GeV0. The system¬ 
atic uncertainties in the effective area of the Fermi-LAT using 
“P7REP_SOURCE_V15” are 10% below 100 MeV, decreas¬ 
ing logarithmically in energy to 5% between 0.316 - 10 GeV 
and increasing logarithmically to 15% at 1 TeV0. We selected 
events with a zenith angle less than 100 ° and that were de¬ 
tected when the rocking angle of the LAT was greater than 
52° to decrease the effects of terrestrial albedo 7 -rays. We 
analysed the 7 -ray data in the direction of MSH 11-6 J A us- 

4 More details found: 

http://www. slac. Stanford, edit/exp/glast/groups/ccmda/la t_Pe rfo rmance. htm 

5 http://fermi.gsfc. nasa. gov/ssc/data/analysis/LAT_caveats. html 
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Figure 1. Left: 2.5° x 1.5° Fermi-LAT count map of the 2 - 200 GeV 7 -ray emission surrounding MSH 11— 61 A (units: counts degree -2 ). The pixel binning is 
0.01° and the maps are smoothed with a Gaussian of width 0.2°. The cyan circles correspond to the background Fermi-LAT sources surrounding MSH 11—61 A. 
The brightest emission seen to the east of MSH 11-61A corresponds to the 7 -ray bright SNR MSH 11-62. The yellow X corresponds to the radio position of 
PSR J1105-6107, while the white X corresponds to the position of ICG J11014-6103 (the lighthouse nebula). One can see that MSH 11-61A is located in a 
complicated region of the 7 -ray emitting sky. 843 MHz MOST radio continuum contours are overlaid in green. Right: 2.5° x 1.5° TS map of MSH 11—61 A. 
The magenta contours corresp onds to the Hi emission of the molecular cloud associated with the remnant as detected by the Southern Galactic Plane Survey 
fMcClure-Griffiths et all20051) . The Hi emission contours range from 66 Kelvin at the edge of the molecular cloud interacting with MSH 11—61A to 100 Kelvin 
at the center of the molecular cloud over seven linearly spaced contour levels. 


ing the Fermi Science Tools v9r33p(Q 

Due to the low count rates of 7 -rays and the large PSF 
of the Fermi-LAT, the standard maximum likelihood fit¬ 
ting technique, gtlike, was used to analyse the 7 -ray emis¬ 
sion of the remnant. Given a specific emission model, 
gtlike determines the best-fit parameters of this model 
by maximising the joint probability of obtaining the ob¬ 
served data given an input model. Gtlike accounts for 
background 7 -ray emission by using diffuse Galactic and 
isotropic emission models described by the mapcube file 
gll_iem_v05_revl. fits and the isotropic spectral 
template iso_source_v05_revl. tx tfl Gamma-ray 
emission from sources found in the Fermi-LAT second source 
catalogue are fixed to their position listed in the catalogue 
and their background contribution is calculated. The Galactic 
diffuse emission arises from the interaction of CRs with the 
interstellar medium and their subsequent decay into 7 -rays, 
while the isotropic component arises from diffuse extragalac- 
tic 7 -rays and residual charged particle emission. 

To improve the angular resolution of the data while 
analysing the spatial properties of the 7 -ray emission of 
MSH 11-6JA, we selected 7 -ray data converted in th e front 
section of the Fermi-LAT with an energy range of 2 - 200 GeV. 
The improvement in spatial resolution in this energy range 
arises from the fact that the 1 -cr containment radius angle for 
front- selected photon events is < 0.3°, while for lower ener¬ 
gies it is much larger. To determine the detection significance, 
position and possible extent of the 7 -ray emission coincident 
with MSH 11-61A we produced test statistic (TS) maps us¬ 
ing gttsmap with an image resolution of 0.05°. The TS is 
defined as 2 log(L ps /L m ,u), where L ps is the likelihood of a 
point source being found at a given position on a spatial grid 
and L nu u is the likelihood of the model without the additional 
source. 

To determine the spectral energy distribution (SED) of the 
7 -ray emission coincident with MSH 11— 6 iA we use events 

6 http://fermi.gsfc.nasa.gov/ssc/data/analysis/software/ 

1 The most up to date Galactic and isotropic emission models can be found 
http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.htmI 


converted in the front section of the LAT that have an en¬ 
ergy of 0.2 - 204.8 GeV. This energy range is chosen to avoid 
the large uncertainties in the Galactic background model that 
arise below 0.2 GeV and to reduce the influence of the rapidly 
changing effective area of the LAT at low energies. We 
model the flux in each of 8 logarithmically spaced energy 
bins and estimate the best-fit parameters of the data using 
gtlike. We also include in the likelihood fit, background 
sources from the 2 4-mo nth Fermi-LAT second source cata¬ 
logue (iNolan et all 120121) that are found within 20° region 
centered on MSH 11— 6 iA. All evident background sources 
were identified in the Fermi-LAT second source catalogue 
and the associated parameters from the catalogue were used. 
We left the normalization of the Galactic diffuse emission, 
isotropic component and the background point sources within 
5 degrees of MSH 11-61A free. For all other background 
point sources with a distance greater than 5 degrees from 
MSH 11-61A their normalisations were frozen to that listed 
in the 2nd Fermi-LAT catalogue. In addition to the statistical 
uncertainties that were obtained from the likelihood analysis, 
systematic uncertainties associated with the Galactic diffuse 
emission were also calculated by artificially altering the nor¬ 
malisation of this background by ± 6 % f rom the best-fit value 
at each energy bin as outlined in lCastro & Slanel (120101) . 

In Figure [T] left panel, we have generated a 2.5° x 1.5° 
count map centered on MSH 11-61A that is smoothed by a 
Gaussian with a width similar to the PSF for the events se¬ 
lected. MSH 11-61A is located in a very complicated re¬ 
gion of the sky. There are many Fermi-LAT sources close 
to the remnant, with 7 -ray bright SNR MSH 11-62 (2FGL 
J1112.1-6040 and 2FGL J1112.5-6105) located ~ 1.2°, and 
PSR J1105-6107 (2FGL J1105.6-6114) located - 0.36°, from 
the remnant. There are four other Fermi-LAT sources in 
the immediate vicinity of MSH 11-6IA (2FGL J1104.7- 
6036, 2FGL J1105.6-6114, 2FGL J1059.3-6118c and 2FGL 
J1056.2-6021), but none of these are coincident with the 
MOST radio contours of the SNR. We define a source re¬ 
gion at the position of MSH 11-6IA to estimate the flux 
from the SNR. Due to the close proximity of PSR J1105-6107 
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Figure 2. The pulse-phase diagram of PSR J1 105-6107 obtained using 
events in an energy range of 0.1 - 300 GeV coming from a 0.5° radius around 
the position of the pulsar. The off-pulse window used for this analysis is de¬ 
fined between 0.00-0.05, 0.22-0.55 and 0.70-1.00 of the pulse phase. Here 
two cycles are shown. 

(J1105) and the low resolution of the Fermi-LAT PSF (~ 1° 
for front events at 68 % containment at ~0.6 GeV), we can’t 
rule out that this pulsar isn’t contributing significantly to the 
observed 7 -ray emission seen in Figure Q] Thus to analyse 
the spatial and spectral characteristics of the 7 -ray emission 
of MSH 11-61A, we have to remove the pulsar contribution. 
As MSH 11-62 is located > 1° away from MSH 11-61 A, it 
is unlikely that it is contributing significantly to the observed 
7 -ray emission. 

2 . 1 . Removing the 7 - ray contribution of PSR J1105-6107 

To avoid contamination from the pulsed emission of PSR 
J1105-6107 (J 1105), we must perform our analysis in the 
off-pulse window of the pulsar light curve. The Fermi-LAT 
collaboration made av ailab le with its second Fermi-LAT cata¬ 
logue of 7 -ray pulsars (lYu et al .120131) the ephemerides of 117 
7 -ray emitting pulsars that had been detected using three years 
of Fermi-LAT data. Due to the continuous observations of the 
Fermi-LAT, lYu et all ( 2013 !) were able to directly determine 
regular times of arrival (TOAs) to produce a precise pulsar 
ephemeris. We use the available ephemeris for J1105 which 
is valid from MJD 54200 (April 2007) to MJD 56397 (April 
2013). As we are interested in analysing over ~70 months of 
Fermi-LAT, we need to check the validity of using the avail¬ 
able J1105 ephemeri s ov er the whole ~70 months. Using the 
Parkes telescope, lYu et al. ( 201 31 analysed J1105 searching 
for timing irregularities over a period of ~ 16 years from Au¬ 
gust 1994 to September 2010 (MJD 49589 - MJD 55461). 
J1105 experienced a glitch at MJD 50417 (November 1996), 
MJD 51598 (February 2000), MJD 54711 (September 2008) 
and MJD 55288 (April 2010). The Fermi-LAT ephemeris of 
J1105 covers the glitches experienced by the pulsar since the 
beginning of the Fermi-LAT mission. As of writing, there 
have been no other reports in the literature that J1105 has un¬ 
dergone glitches since April 2010. We also tested whether the 
pulsar is noisy in 7 -rays, as this can also indicate irregulari¬ 
ties in the pulsar rotation that have not been incorporated in 
the ephemeris. We produced a pulse phase vs. events plot 
during the period of validity of the ephemeris (MJD 54200 
to MJD 56397) and for the period between the end time of 
the ephemeris and the end date of our data (MJD 56397 to 
MJD 56824) to see if the pulse peaks that we observe in Fig¬ 


ure|2]disappear due to noise of the pulsar. We find that during 
both time periods we obtain the same pulse phase profile. As 
no glitches have been detected as of writing and the pulsar 
has not been noisy beyond MJD 56397, we assume that this 
ephemeris is valid for our whole data set. 

The 7 -ray photons were phase-folded using this ephemeris 
and the pulse phases were assigned to the Fermi-LAT data us¬ 
ing the Fermi-LAT TEMP02 plugin provided by the Fermi- 
LAT collaboratiorfl. This plugin calculates the rotational 
phase of the pulsar for each photon arrival time in the Fermi- 
LAT data using the barycentric dates of each event. Using 
ft select, we remove the pulse and use only the 7 -ray photons 
in the off-pulse window (defined by the 0.00-0.05, 0.22-0.55 
and 0.70-1.00 pulse phase intervals) to perform our spectral 
and morphological analysis. In Figure [2] we have plotted the 
pulse-phase diagram of J1105 obtained in the 0.10-300 GeV 
energy range using 0.5° radius around the position of J1105. 

2 .2. TS Map 

In Figure [j] right panel we present a TS map of 
MSH 11-6JA using the off-pulse 7 -ray data. This was cal¬ 
culated using gttsmap over an energy of 0.2 - 2.0 GeV and 
using front events only. In addition to the diffuse Galac¬ 
tic background components, we include in the background 
model the Fermi-LAT sources associated with MSH 11-62 
(2FGL J1112.1-6040 and 2FGL J1112.5-6105) and the four 
sources in the immediate vicinity of the remnant (2FGL 
J1104.7-6036, 2FGL J1105.6-6114, 2FGL J1059.3-6118c 
and J1056.2-6021). The TS map suggests that there is signif¬ 
icant 7 -ray emission coincident with MSH 11-61A and the 
MC associated with the remnant, as highlighted by the ma¬ 
genta contours. The peak of the 7 -ray emission is found at 
a best fit position of (a y 2 ooo, ^ 2000 ) = (1 l /, 01 m 29',-60°5529), 
placing it outside the SNR boundary, but consistent with be¬ 
ing located along or inside the western limb given the angular 
resolution of the Fermi-LAT. The emission is detected with 
a significance of ~ 5<r, with the 7 -ray emission at the center 
of the remnant producing a significance of ~ 4cr. One can 
see in Figure Q] right panel that the contribution of MSH 11- 
62 and the other Fermi-LAT sources surrounding the remnant 
have been modelled out, while the contribution from the pul¬ 
sar (J 1105) has been gated out successfully. 

2.3. 7 -ray spectrum 

The 7 -ray spectrum of MSH 11-6 JA is shown in Figured 
with the statistical errors plotted in black and systematic er¬ 
rors plotted in red. For energies above 6.40 GeV, only flux 
upper limits have been determined and are plotted as blue tri¬ 
angles. Additionally, the best fit power law and exponential 
cut-off power law models are plotted as the green dotted line 
and purple dot-dashed line respectively. The 7 -ray spectrum 
can be fit using a simple power law with a spectral index of 
2.75^o{}g, giving a reduced x 2 ~ U An exponential cut-off 
power law with E cut = 4.20^Qgg GeV and spectral index of 
2'49 ^q [9 can fit the spectrum equally well giving a reduced 
X 2 ~ 0.8 for the fit. For an energy range of 0.1-100 GeV, the 
integrated flux is (4.2^ ^) x 10 ~ u erg cm - 2 s -1 , assuming the 
power law fit. Using a distance of 7 kpc, the luminosity of this 
7 -ray source in this energy range is ( 2 . x 1035 erg s '. 

3. SUZAKU OBSERVATIONS OF MSH 11-6 JA 

8 http://fenni.gsfc.nasa.gov/ssc/data/analysis/user/Fenni_plug_doc.pdf 
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Figure 3. The Fermi—LAT 7 -ray spectrum of MSH 11— 61 A. Statistical er¬ 
rors are shown as black error bars, systematic errors are plotted as red error 
bars, while the upper limits are plotted as blue triangles. The simple power 
law model with T = 2.75 ^‘qI and exponential cut-off power law model with 
T = 2.49 ^']q and E C ut = 4.20^‘gJ GeV are shown as the green dotted and 
magenta dot-dashed line respectively. 


MSH 11-61A was obse rved with Suzaku usin g X-ray imag¬ 
ing spectrometers (XIS) iKovam a et al.ll2007l) on 2011 July 
25th for ~ lllks (ObsID 506061010). For this observation 
only XIS(fl XIS1 and XIS3 observations are available as 
XIS2 has not been functi o nal si nce November 200(0 Re¬ 
cently, iKamitsukasa et al] (120151) presented their analysis of 
this Suzaku observation. They found recombining plasma in 
the center and in the northwest of the remnant which has en¬ 
hanced abundances and a temperature of ~0.5 keV, while ev¬ 
erywhere else the X-ray emission arises from an ionising ISM 
component with a temperature of ~0.6 keV. In section 5.2.1, 
we estimate the density of the 7 -ray emitting material based 
on our Fermi-LAT spectrum (Figure [3}- To test whether this 
inferred density agrees with other observations, we have re¬ 
analysed the Suzaku data in order estimate the density of the 
surrounding environment. 

For our analysis we used the standard tools of HEASOFT 
version 6.16. We reprocessed the unfiltered public data us¬ 
ing aepipeline (version 1 . 1 . 0 ) and use the current calibra¬ 
tion database (CALDB) available as of 2014 July 1st (version 
20140701). Following the standard screening criterion, we 
filtered hot and flickering pixels, time intervals corresponding 
to Suzaku passing the South Atlantic Anomaly and night-earth 
and day-earth elevation angles less than 5° and 20°, respec¬ 
tively. We utilised events that had a grade of 0, 2, 3, 4 and 6 
only. The total exposure of our observation is 111 ks for each 
of the XIS detectors. We extracted the spectra and images of 
the remnant from the 5x5 and 3x3 editing mode event files 
using XSELECT version 2.4. For the spectral analysis we 
generated the redistribution matrix file (RMF) and ancillary 
response files (ARF) using XISRMFGEN and XISARFGEN 
respectively. To analyse the spectral data we used the X- 
ray spectral fittin g package (XSPEC) version 12.8.2 q with 
AtomDB 3.0. lPl (iSmith et ahfoOQlHFoster et al. 1120121) . 

9 It is also important to note that a fraction of XIS0 has 
not been functional since 2009 June 23rd due to the dam¬ 
age caused by a micro meteorite. For more information see: 
http://www.astm.isas.jaxa.jp/suzaku/doc/suzakumemo/suzakumemo-2010- 
01. pdf. 

10 http://www.astro.isas.jaxa.jp/suzaku/doc/suzakumemo/suzakumemo- 
2007-08.pdf 

11 http://heasarc.nasa.gov/docs/suzaku/pmcessing/criteria xis.html 

12 AtomDB 3.0.1 can be downloaded here: 
http://www. atomdb. org/download.php 


Table 1 

Best tit parameters for our background model. 


Component 

Parameter 

Value 

Cosmic X-ray background 

N//(x 10 22 cm' 2 ) 

r 

1 47+i.u/ 

1.40 (frozen) 

Galactic ridge emission 

N//(x 10 22 cm' 2 ) 
kT (keV) 
Abundances 

0 61 mus 
uol -0.06 

0.22 ± 0.02 
0.20 (frozen) 

Galactic halo 

N//(x 10 22 cm' 2 ) 

kT (keV) 

Ne 

Mg 

Si 

1 64 +tuJ 

*■ -0.10 

0.58 ±0.02 

9 71 +O. 6 O 

-0.18 

9 48+0.34 


Reduced (dof) 

1.80(1602) 


Note: All uncertainties correspond to lcr errors. 

3.1. Spectral analysis of the individual annulus and 
rectangular regions 

We extracted spectra from a central circular region defined 
by region 1 in Figure 4 (right panel) and 6 annular regions of 
width 0.82' to cover the central X-ray emission of the rem¬ 
nant (regions 2 - 7 in Figure 4, right panel). The radial size 
of these regions was chosen to be the same size as the angu¬ 
lar resolution of Suzaku (~ 0.8'). These regions were chosen 
to fully enclose the bright central X-ray emission of the rem¬ 
nant, which is quite symmetric in nature. We also extracted 
spectra from three rectangular regions (regions 8 - 10 in Fig¬ 
ure 4, right panel) that are not covered by the annulus regions, 
to enclose protrusions in the northwest, southeast, and east. 
Annular regions were chosen in order to characterise radial 
variations in the brightness, temperature, ionisation state and 
elemental abundances of the remnant, all of which are im¬ 
portant for understanding the nature of the mixed morphol¬ 
ogy. Our choice of reg ions differs from those chosen by 
IKamitsukasa et al] (120151) . who also analysed the Suzaku ob¬ 
servation of MSH 11-61A. They extracted spectra from five 
regions that do not enclose the full X-ray emission from the 
remnant - a central region that corresponds to our regions 1 , 
2, and 3; a northwest region that encompasses our region 8 
and a northwestern portion of region 7; and NE, SE, and SW 
regions that cover sectors of our annular regions and also en¬ 
compass our region 9 in the east. All spectra were grouped 
by 20 counts using the FTOOLS command grppha and all 
available XIS detectors were used. 

To estimate the background, we extract data from the full 
field of view of the XIS of our observation, excluding the cali¬ 
bration regions and the emission from the remnant. The back¬ 
ground spectrum consists of two major components, the non 
X-ray background and the astrophysical background which is 
made up of the cosmic X-ray background, the Galactic ridge 
X-ray emis sion and the Galactic halo. We use xisnxbgen 
dTawa et all 12008b to generate a model for the NXB which 
we then subtract from our background spectrum. Similarly, 
we subtract a model for the NXB from our spectra obtained 
fr om the regions shown in Figure |4] right panel. Similar 
to IKamitsukasa et all (120151) . we model this NXB subtracted 
spectrum. We fix the cosmic X-ray background power law 
component to that of iKushino et all (2002 ). and use a sin¬ 
gle apec model with a temperature and surface brightness 
(~ 1.0 x 10 -12 erg cm ' 2 s ' 1 deg ~ 2 for 0.5-2. 0 keV) similar 
to that obtained by iHenlev & Shelton! (120131) to define our 
Galactic halo component. We use a single low tempera¬ 
ture apec model with subsolar abundances frozen to that of 
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Figure 4. Suzaku XISO image of MSH 11—6iA in the 0.5-7.0 KeV energy band. The image has been smoothed with a Gaussian function of width 0.1’ and has 
a logarithmic scaling applied to it. Left : The green contours correspond to the 843 MHz MOST radio continuum emission used in Figure Q] Right : Overlaid in 
green are the regions we use for spectral extraction, described in Section 3.1. 


Table 2 

Spectral fits for all 10 individual regions. 


N h (10-- cm--) 


Te r (10 12 s cm 4 ) Reduced \ 2 


Region 


kT (keV) 

0 43+<>-U3 

-QQ 2 
n 42+0-02 
u -QQ 2 

0 4Q+0-03 

u.^u_oo 2 

n 34+0.02 
^ -0.09 

0 27 +0 06 
u , z / - 0.01 
0 30 + 0-03 

U. JU _0 02 

0.363-2? 


kT iml (keV) 


Ne 


Mg 

90+0.52 

q 4 + 8 :?§ 
■ -0.26 
qo+0.33 
•^ J -0.28 
4Q+0.33 
-0.28 
67 +0.51 
• o / - 0.20 
+1.07 

59 + 8:12 

,OJ_I 02 

69+1 ■ 67 
■ ov -l.n3 


Si 

+1.07 


s 

i+0.!14 


' s cm 

1 59+U-19 

1.29 

1 26+™ 

1 34+0.04 

l-283?f 

I no+0 08 
+■82 
0 9O+-05 
+ Q 6 
0 + 2 + 0-22 
u ~ z - 0.11 
o 19+0.39 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


,t 8 :B 

,78:11 

1:11 


1.76; 
1.20 
1.32 
1.81 

i 87+»:» 

1 fifi+O - 04 

1 -0.08 

1 67+ 0 - 21 
227^:15 
20078:11 

-0.30 

2.46+0-38 


0.42 
0.80_0 27 
0.73+O- 24 


;t8:S4 


5 

5 

5 

5 

5 

5 

5 

5 

0.0808 

0.0808 


0 23 +0 - 12 
o 35 + 8 :?? 

U.JJ_0 1 2 


.11 


6.32 

4.96 

4.65 




1.92 

2.51 

2.62 


'78:18 


5 9 ]+0.74 
J Z -0.75 

3 69+0- 67 

J O -0.41 

7 05+*- 92 

, u -1.30 

4 87+ LU 

°- yi -1.96 

6.67+f?? 


0.11 

0.14 


78:?? 


5.69^.73 3.72^- 4 J 

\52 


3 34+1.52 
9 i 9+0.46 
z - az -0.44 

4 97+1-42 

991^:11 

306f,l 

2 . 867‘- 4 3 


,78:11 


0.79 

0.94 

1.02 

1.08 

1.06 

1.10 

1.14 

0.92 

0.84 

0.94 


Note: All uncertainties correspond to the 90% confidence level. 


iKaneda et al J (1 1 9971) to define the Galactic ridge emission. W e 
also use the Wilms et al. abundance table dWilms et al.ll2000|) . 
Our best-fit parameters for our background and their uncer¬ 
tainties are given in Table [j] 

To model the X-ray emission from the remnant we used 
a non-equilibrium ionisation (NEI) collisional plasma model, 
VVRNEI, which is characterised by a final (kT) and initial 
electron temperature (kT init ), elemental abundances and a sin¬ 
gle ionisation timescale (r = n e t). This allows one to model 
a plasma that is ionising up to collisional equilibrium from 
a very low initial temperature kTj n i t , mimicking the standard 
NEI/VNEI model that is commonly used in the literature. 
Additionally, the RNEI/VVRNEI model can reproduce a re¬ 
combining (overionised) plasma where one assumes that the 
plasma starts in collisional equilibrium with an initial temper¬ 
ature kTi n n that suddenly drops to its final temperature kT. For 
our analysis, the column density, ionisation timescale, nor¬ 
malisation, and final temperature were left as free parameters. 
Due to the strong emission lines from Mg, Si, and S the abun¬ 
dances of these elements were also left free. Additionally we 


also let Ne and Fe be free parameters for regions 2 - 3, as we 
found that varying these significantly improved the fit. All 
other elemental abundanc es w ere frozen to the solar values 
reported by IWilms et al.l (120001) . The foregrou nd absorbing 
column density Nh was modelled using TBABS (1 Wilms et al.l 
120001) . Figure [7] shows an example of the X-ray spectrum of 
MSH 11-6J A as extracted from region 3. The spectra derived 
for each region shown in Figure Q] right panel all have similar 
features to the spectrum shown in Figure|7] 

We found that for regions 1-8 the fit favoured an initial 
temperature larger than the final temperature, while regions 
9 and 10 favoured an initial temperature smaller than the fi¬ 
nal temperature. When left free, these initial temperatures 
would hit the upper (or lower) limits of this parameter and 
the associated abundances we obtained for our fits were un¬ 
realistically high (abundances of [Mg], [Si] > 10 relative to 
IWilms et al.ll2000l) . The ionisation timescale for all regions 
was r ~ 10 12 s cm -3 . 

To investigate the sensitivity of our fits to values of kT m i t , 
we simulated spectra with similar counting statistics to those 
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Figure 5. The sensitivity of our fits to changes in based on Ax 2 - Here 
the red plots correspond to kTj n it = 5 keV, and the blue corresponds to kT m i, 
= 20 keV. The solid, dotted and dot-dashed lines correspond to r = 10 10 s 
cm -3 , 10 11 s cm -3 and 10 12 s cm -3 respectively. We assume kT = 0.5 keV 
and riH = 10 22 cm -2 for all spectra. From this plot we derive an upper limit 
for kTi n i t of 2 - 5 keV corresponding to a shock velocity of 1300 - 2100 km/s, 
for which we then use in our fits. 


from our regions of investigation, using kT = 0.5 keV, Nh = 
10 22 cm -2 , and kT init - 5 keV or 20 keV. We considered 
cases for logr = 10, 11, and 12. We fit each spectrum to 
a TBABS*VVRNEI model and then investigated the effect 
of freezing kT lni , values over a range from 0.1 to 100 keV. 
The results are illustrated in Figure [5] where we plot Ay 2 
vs. kTi n i t for spectra with actual kT lnn values of 5 keV (red) 
and 20 keV (blue). Here the solid, dotted, and dashed lines 
correspond, respectively, to logr = 10, 11, and 12. For low 
ionisation timescales the resulting fits are quite sensitive to 
the kTi n i t value, while for timescales similar to that found in 
MSH 11-61A (r « 10 12 s cm -3 ), the fits are insensitive to 
kti n it (i.e. the Ay; 2 vs. k7}„„ plot plateaus) for temperatures 
greater than 2-5 keV. This is because above ~ 2 - 5 keV, it 
is only the emission from Fe and Ni that is significantly im¬ 
pacted by higher kTnit values, because all other abundant ions 
are fully stripped at these temperatures. Our observations do 
not have enough counts around the Fe and Ni lines to pro¬ 
vide sensitivity to such an effect (though longer observations, 
particularly with higher resolution, would provide such sensi¬ 
tivity; see below). 

Based on where Ay 2 vs. kT init plot plateaus for r ss 
10 12 s cm -3 , we fix value of kT in j, at 5 keV for all regions 
other than 9 and 10. For the latter regions, whose fits indicate 
kTjnit < kT), we fix the initial temperature at the minimum 
available value for the VVRNEI model (80.8 eV). In Table 
[2] we list the best fit parameters for each of our spectra. All 
uncertainties correspond to the 90% confidence level. 

A recombining plasma that has an initial temperature of 
2-5 keV implies that the shock must have had a veloc¬ 
ity of ~ 1300 - 2100 km s -1 , assuming electron-ion equilib¬ 
rium. This high initial velocity suggests that the recombin¬ 
ing plasma could be the result of the SNR shock front ini¬ 
tially expanding into a dense CSM, reaching a high ionisa¬ 
tion state corresponding to the high initial expansion velocity 
(and, thus, shock temperature). Such a scenario might result 
from expansion into an r -2 density profile characteristic of a 
stellar wind, with subsequent expansion into the lower den¬ 



Energy (keV) 


Figure 6. A 20 ks background subtracted spectrum simulated for ASTRO-H 
for a recombining plasma that has an initial temperature of 2 keV (black) and 
one that has an initial temperature of 50 keV (red). The model was based on 
that obtained for region 8 but with the initial temperature set to 2 keV and 50 
keV. With this observation one can easily differentiate between a recombining 
plasma that has two different initial temperatures. 


sity regi ons resulting in rapid cooling, leav ing an overionised 
plasma (lltoh & Masailll989fclMorivall2012l) . 

Even though our current data are unable to differentiate be¬ 
tween a plasma that has an initial temperature of 2 - 5 keV 
or one that has initial temperature greater than this, with the 
launch of ASTRO-H we will be able to differentiate between 
these two cases. In Figure [6] we have plotted a simulated 
background-subtracted spectrum that we would obtain in a 
20 ks observation using the calorimeter on ASTRO-H for a 
plasma that has an initial temperature of 2 keV and one that 
has an initial temperature of 50 keV, assuming our fit param¬ 
eters for region 8. The background spectrum for ASTRO-H 
was obtained from SIMX simulation^. Here the black spec¬ 
trum corresponds to a plasma that has an initial temperature 
of 2 keV (or an initial velocity of ~ 1300 km s -1 ), while the 
red spectrum corresponds to a plasma that has an initial tem¬ 
perature of 50 keV (or an initial velocity of ~7000 km/s, typ¬ 
ical of the high initial expansion speed of an SNR). One can 
see that with a 20 ks ASTRO-H observation, we could easily 
differentiate between two plasmas that have different initial 
temperatures. 

Our spectra from MSH 11 -6 / A are best described by 
a recombining plasma model with the exception of emis¬ 
sion from regions 9 and 10, in the eastern and southeast¬ 
ern outskirts, where an ionising plasma is observed. Our 
re sults for the central region s (1-3) agree well with those 
of iKamitsukasaetahl d2015l) , who also find a recombining 
plasma for their southwestern region, in agreement with our 
results . In contrast, for their sout heastern and northeastern re¬ 
gions, IKamitsukasa et al. ( 201 5) obtain best fits for an ionis¬ 
ing plasma. Given that these regions combine emission from 
the eastern and southeastern protrusions (our regions 9 and 

13 https://hea-www. hcin’ard. edu/simx/ 
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1 2 

Energy (keV) 


Figure 7. The Suzuku XISO, XIS1 and XIS3 spectrum extracted from region 3 in Figure[4] These spectra were fitted simultaneously using our background model 
described in Tableland an absorbed VVRNEI model as described by the parameters listed in Table[2] The spectra are overlaid with the best fit model with their 
chi 2 residuals. The spectra derived for each region shown in Figure fright panel all have similar features as the spectrum shown here. 


10 ), for which we also observe an ionising plasma, with emis¬ 
sion from the outer symmetric portion of the SNR, where 
we observe a recombining plasma, the combination of two 
components may explain the partial discrepancy. Conversely, 
iKamitsukasa et all d2015b report an ionising plasma for their 
southwestern region, which covers regions for which we ob¬ 
tain a recombining plasma. This may result from our annular 
regions averaging over multiple components. 

All regions have an ionisation timescale of ~ 1.0 x 
10 12 cm -3 s, indicating that the X-ray emitting plasma 
across the whole remn ant is close to ionisation equilibrium 
(ISmith & Hughesf2010b . The temperature of the recombining 
components ranges from 0.27 keV to 0.43 keV near the bright 
central X-ray emission. Interestingly, the regions that are best 
described by an ionising plasma have the highest temperatures 
out of all the regions that we analysed, with temperatures of 
0.80 keV and 0.73 keV respectively. On average, _our derived 
tempe r atures are lower than that derived bv lKamitsukasa etafl 
(120 1 51). iG arcfa et all (120 1 2b in their XMM—Newton analysis, 
and lSlane et all (120021) in their ASCA analysis. 

We find in all regions strong emission lines coming from 
Mg, Si and S and all regions require super-solar abundances 
of these elements. The enhancement of elemental abun¬ 
dances is observed in many MM SNRs and indicates that 
some of the X-ray emission we are observing arises from 
ejecta that have been dispersed throughout the remnant and 
been mixed with t h e swe pt up shocked material. S imilar to 
IKamitsukasa et all (120151) . and IGarcfa et al.l (1201 2b we also 
find an underabundance of Ne and Fe in regions 2-3. 
Unlike, IKamitsukasa et all (120151) we do not find evidence 
for overabundance of Ar or the underabundance of O sug¬ 


gested by IGarcfa et al.l (12012b . When we freed these param¬ 
eters we found that they do not significantly improve our 
fit, thus we kept them at solar. Our estimated abundances 
are sli g htly higher t h an tha t der ived bv IKa mitsukas a et~afl 
(1201 5b . IGarcfa et al.l (1201 2b and ISlane et all (120021k This 
discrepancy arises from the fact_that_in our analysis we 
use the abundance table by I Wilms et all (12000b and the 
newly upda ted ATOMDB 3.0. 1, while IKamitsukasa et all 
(1201 5b and IGarcfa et al.l (1201 2b use the table derived by 
I Anders & Grevessel (1 1989b and ATOMDB 3.0 and ATOMDB 
2 .0.2 respectively. 

Our derived column density towards MSH 11-61A ranges 
between (1.20j{5'[ 2 -2. 46 ^ 037 ) x 10 22 cm -2 . The column den¬ 
sity is highest in regions 8-10 which directly interacts 
with the surrounding environment. Our e stimates for Nh are 
higher than the colu mn density derived bylSlane et al.l (12002b . 
IGarcfa et al.l (1201 2b and IKamitsukasa et al.l (1201 5b . This dis¬ 
crepancy arises from the fact that we use a different abun¬ 
dance table and a different absorption model. 

3.2. Deriving the density of the X-ray emitting gas 

The density of the X-ray emitting gas was calculated from 
the normalisation of the VVRNEI models using n e = 1.2 nn- 
We estimate the volume for each region by taking an area 
equivalent to the extracted SNR regions shown in Figure [4] 
right panel and projecting this area through a filled sphere. 
The estimated density n « 1.1 nn (assuming ISM abundances) 
is listed for each region in Tabled 

The inferred density ranges from n = (0.66^q 23 — ^ .87l!l| ||) 
d-o.5f-o.5 cm -3 anc j j s highest in region 8 which is coincident 
with the location of the dense molecular cloud found towards 
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Table 3 

Number density of the X-ray emitting material estimated from the best fits 
of the 10 different regions shown in Figure[4] 


Region 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


(d-0.5 fD.5 cm 


o 66+h-H 

u.oo_ 0 23 

0.85!; 

1.28 + ltl 


2.43! 

2 . 68 ! 


;l:3I 
-1.20 
0 ai+1.28 

i -^ 1 -1.15 

5 . 87 !^? 

1 77+^32 
1,1 1 -0.93 

±£±zUi 


the west. The density is lowest towards the center of the rem¬ 
nant where the brightest X-ray emission is located, while the 
eastern part of the remnant has a density that is intermediate 
of the central regions of MSH 11-61A. Our density estimates 
for the bulk of the remnant are consistent with the densities 
derived bv ISlane et al.l (12003 ) who attempted to reproduce the 
observed temperature and brightness_profiles of the remnant 
using the cloudy ISM model by I White & Lond (119911) and a 
hydrodynamical model that traces the evolution of the rem¬ 
nant, while incorporating the effects of thermal conduction. 

4. THE ORIGIN OF THE THERMAL X-RAY EMISSION 

The total X-ray emitting mass in MSH 11 -6 / A is given 
by M x = lAn H m H fV, where m H is the mass of the hydrogen 
atom, V is the volume from which the emission is observed, 
and / is the filling factor. Using the estimated volumes and 
derived densities for regions 1-7, we sum the masses to obtain 
Mx ~ 480 d '!'' 2 !' 1 ! 2 M-... This is co mparable to the swept up 
mass derived bv ISlane et al.l (120021) . 

The enhancement of Mg, Si and S abundances throughout 
the remnant suggests that the observed X-ray emission arises 
in part from supernova ejecta. Assuming that all ejecta have 
been shocked, we can estimate the mass of the ejecta com¬ 
ponents based upon the measured abundances: M, = [(a, - 
1 )/ 1 .4 ](ni/n H )(mi/m H )Mx where M t is the ejecta mass of 
species i, a, is its abundance relative to ISM abundances, 
as listed in Table 2, m, is the atomic mass, and n,/n h is 
its ISM abundance relative to hydrogen. We find that, us¬ 
ing the average of the measured abundances, the total ejecta 
masses of Mg, Si, and S are, respectively, 0.37^ 5 / 2 / 1 ! 2 Mq, 
0.80and 0.26d 5 ' 2 f l / 2 M Q . However, we note 
that the abundances for Ne and Fe are both lower than ISM 
values, meaning that we have no evidence for ejecta compo¬ 
nents for these species, and suggesting caution in interpreting 
all of the abundances. Taken at face value, however, the Mg, 
Si, and S ejecta m ass estimates are consiste nt with a progeni¬ 
tor mass > 25Mq ([Thielemann et al.lll996l) . 

Recombining plasma can arise from two main scenarios: 
thermal conduction which is the rapid cooling of electrons 
due to the interaction of the ho t ejecta w ith the cold, dense 
surrounding environment (ICox et all 11999 ): or adiabatic ex¬ 
pansion which can occur when the SNR shockwave expands 
through a dens e circ umstellar medium into a low density ISM 
(lltoh & MasailU989l) . 

For thermal conduction to be the more likely scenario, the 
recombining plasma is expected to be coincident with the 


location of the molecular cloud, there should be a tempera¬ 
ture decrease towards the molecular cloud and one would ex¬ 
pect the thermal conduction timescale t co „d to be less than or 
comparable to the age of the remnant. We find recombin¬ 
ing plasma in regions that are directly interacting with the 
molecular cloud (Regions 7 and 8 ), and we do see a slight 
temperature decrease towards the molecular cloud based on 
the annulus regions. The thermal conduction timescale is 
given by dSpitzeril 1962t [Zhou et al.ll2014l) t con j ~ kn e lj/K ~ 
56(n e / \cm~ 2 )(lj / l0pc) 2 (kT /0.6keV) -5 / 2 (Z«A/32)kyr, where 
n e is the electron density and is calculated from our best-fits 
listed in Table [2] lj is the scale length of the temperature gra¬ 
dient, k is Boltzmann’s constant, k is the thermal conductivity 
for a hydrogen plasma and In A is the Coulomb logarithm. 
Assuming a distance of 7 kpc to the remnant and a radius 
of ~ 5.11', we calculate the length of the temperature gradi¬ 
ent to be ~ 3.6 x 10 19 cm. Using the temperature and density 
of region 7 (see Table [7] and Table [3 respectively), the ther¬ 
mal conduction timescale is estimated to be ~ 334 kyr. This 
is ~H 6 times greater than the age of the remnant (~ 10-20 
kyr), making it unlikely that the overionised plasma arises via 
thermal conduction. 

Another possibility is that the recombining plasma arises 
from adiabatic cooling. To calculate t recom b, we use the best 
fit ionisation timescale for region 1-8 listed in Table [2] and 
divide these by the electron density of each region. We obtain 
a recombining timescale between ~ (2— 125) d^' 5 f “°' 5 kyr, 
which is comparable to the age of MSH 11-61A, making this 
scenario t he most likely. This is con sistent with the results re¬ 
ported by iKamitsukasa et al . (201 5) and the velocity implied 
by our upper limit for kT init . 

5. THE NATURE OF 7 -RAY EMISSION 

5.1. Pulsar contribution and the Integral source ICG 
J11014-6103 

Pulsars that are detected within the Fermi-LAT energ y b and 
(see the second Fermi-LAT Pulsar catalogue by lAbdo et all 
120131) . have spectra that is well characterised by a power law 
with an exponential cut-off of 1 - 5 GeV. As the 7 -ray spec¬ 
trum of MSH 11-61A can be described using an exponential 
cut-off of E cu , ~ 4.2GeV, we still need to consider the sce¬ 
nario that the 7 -ray emission we observe arises from a nearby 
pulsar other than J1105-6107. 

Using the Australian Tel escope Natio nal Facility (ATNF) 
Pulsar Catalogue dManchester et al.ll2005l) . there are 9 pulsars 
including J1105-6107 within 5°, whose spin down power is 
sufficient to produce the 7 -ray flux of MSH 11-61A. All of 
these pulsars, except for J1105-6107, are >1° from the cen¬ 
troid of the 7 -ray emission making it unlikely that any of these 
pulsars are contributing significantly to the observed emission 
of MSH 11-6IA. As we removed the contribution of J1105- 
6107 from the 7 -ray data as described in Section 2.1, we can 
also rule ou t its c ont ribu tion. 

Recently, iPavan et ail (120141) investigated the nature of the 
X-ray and radio emission of the INTEGRAL source ICG 
J11014-6103, which they call the lighthouse nebula. In X- 
rays this nebula exhibits a prominent jet-like feature that 
is perpendicular to an elongated cometary tail, and a point 
source. The source of this X-ray structure is a neutron star 
travelling supersonically and we have plotted its position as 
the white X shown in Figure Q] This neutron star has a spin 
down power of /! ~ 10 37 erg s -1 . As the 7 -ray luminos¬ 
ity of MSH 11-6IA is 2.5 x 10 35 erg s“‘, the PSR of ICG 
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J11014-6103 would require an efficiency of L 7 /E ~ 2.5% to 
produce the observed 7 -rays, which is plausible. In an at¬ 
tempt to disentangle the likely source of the 7 -ray emission, 
we have plotted as the magenta contours in Figure Q] right 
panel, the Hi contours of the molecular cloud associated with 
MSH 11-61 A. If one assumes that the pulsar of ICG J11014- 
6103 can produce significant 7 -ray emission, we would ex¬ 
pect the detection significance peak to be skewed towards the 
position of ICG J11014-6103 instead of in the direction of the 
remnant and MC as is observed. Thus even though we can¬ 
not rule out that ICG J11014-6103 is contributing to the 7 -ray 
emission in Figure [3] the association of the molecular cloud 
and the detection significance in this region suggests that the 
emission most likely arises from the interaction of the SNR 
with the molecular cloud, rather than ICG J11014-6103. 


5.2. Modelling the broadband emission of MSH 11—61A 

To investigate the nature of the broadband emission from 
MSH 11-61A we use a model that calculates the non-thermal 
emission produced by a distribution of electrons and protons. 
The 7 r° de cay m odel is based on the proton-proton interaction 
model by iKamae et all (120061) . with a scaling factor of 1.8 5 
for helium and heavy nuclei as suggested by iMor j| <m§. 
Th e leptonic emission mod els are based on tho se presented 
by iBaring et~akl (119991) and iBvkov et~afl (120001) for the syn- 
chrotron/IC and non-thermal bremsstrahlung emission mech¬ 
anisms. We assume a spectral distribution of our accelerated 
particles to be: 


aly i -a, 

=a,p 'exp 
dp 



( 1 ) 


where i is the particle species, a, is the normalisation of the 
particle distribution, a, is the particle distribution index which 
is equal to (l-r)/2, where I’ is the photon index and /?(>,■ is 
the exponential particle momentum cut-off. This distribution 
is transformed from momentum space to energy space such 
that the exponential cut-off is defined by an energy input, Eqi- 
The sum of the integrals of each spectral distribution is set to 
equal the total energy in accelerated particles within the SNR 
shell, Ecr = cEsnr , where e is the efficiency of the SNR in 
depositing energy into cosmic rays. 


5.2.1. Hadronic origin of the observed y-rays 

In Figure [ 8 ] we have plotted the model fits to the broad¬ 
band emission of MSH 11-61A. The radi o spec tru m is a 
combination of mu l tiple o bserv ations bvlMilne et al.l (1 1 9891) . 
IWhiteoak & Greenl (119961) . and iFilipovic et al.l (120051) . The 
X-ray upper limit was derived by fitting a power law 
with a photon index similar to that of RX J1713.7-3946 
(lUchivama et al.l [20031) . Kepler and RCW 86 (iBamba et al.l 
120051) (T = 2.3), to our models of the Suzaku data. The 
upper limit corresponds to the flux in which the additional 
non-thermal component begins to affect our reduced yf ■ 
The solid magenta line corresponds to the 7 r°-decay model 
that adequately reproduces the observed 7 -ray spectrum of 
MSH 11-61 A. We have also plotted as the purple dot-dashed 
line the synchrotron model that sufficiently reproduces the ra¬ 
dio spectrum, assuming an electron-proton ratio ( K ep ) of 0 . 01 , 
while the IC model falls below the plotted axis. For complete¬ 
ness we have also plotted the non-thermal bremsstrahlung 
contribution as the orange dashed line. Table [4] lists the pa¬ 
rameters which reproduce the 7r°-decay, synchrotron emis¬ 


sion, IC and non-thermal bremsstrahlung models plotted in 
Figure [ 8 ] 

A 7 r°-decay model arising from a proton distribution with 
a power law index of a p = 4.34 and a cut-off energy of 
E{’ ~ 6.05 GeV, adequately described the 7 -ray spectrum of 
MSH 11-61A. The cut-off energy of the proton spectrum de¬ 
rived in this model is much smaller than the Te V cut-off one 
would expect for protons ( Reynolds 2008 ). This could indi¬ 
cate that due to the high density of the surrounding environ¬ 
ment, efficient CR acceleration is suppressed allowing accel¬ 
erated parti cles to escape the emission volume dMalkov et al.l 

l2mHl20l2h . 

As non-thermal X-ray emission has not been observed from 
MSH 11-61A, we are unable to constrain the cut off energy 
of the electron population. Thus to model the radio emission 
of MSH 11-61A we assume the electron distribution has the 
same cut-off energy as the proton distribution. We are able 
to reproduce the radio spectrum using an electron distribution 
that has a power law index of a e = 3.15 and a magnetic field 
of 28/iG. The magnetic field implied by the synchrotron mod¬ 
elling is larger than the magnetic field of the ISM (~ 3-5/iG). 
This enhancement could arise from magnetic field amplifica¬ 
tion due to the compression of the surrounding medium by the 
SNR shock-wave. 

In the 7 r°-decay model, the calculated 7 -ray flux is pro¬ 
portional to the ambient density of the surrounding ambient 
medium and the total proton energy. Assuming a conserva¬ 
tive upper limit of 40% of the total supernova explosion en¬ 
ergy goes into accelerating CRs, we can estimate the density 
of the 7 -ray emitting material. For our 7 r°-decay model of 
MSH 11-61A we obtain an ambient density of 9.15 cm -3 . 
Similar to many other SNRs that exhibit hadroni c emi ssion 
(e.g W41, MSH 17-39, G 337.7-0.1: ICairo et al.llExinl) : Kes 
79: lAuchettl et al.l (120141) ). this density is much larger than 
the ambient density estimate derived from our X-ray analysis 
(see Table [3j. This discrepancy could arise from the SNR 
shock-wave interacting with dense clumps of material that 
are c old enough such that do not radiate sig nificantly in X- 
rays (ICastro & Sland2010Ulnoue et al .l20~i~2lr If these clumps 
have a high filling factor, then the densities that we derive in 
our X-ray analysis would underestimate the mean local den¬ 
sity. Our inferred ambient densities as well as the association 
of MSH 11-61A with a molecular cloud towards the west of 
the remnant supports the conclusion that MSH 11-61A is in¬ 
teracting with dense material that does not radiate in X-rays. 

An alternative scenario is that the enhanced 7 -ray emission 
arises from highly energetic particles escaping the accelera¬ 
tion region and are interacting wi t h den s e gas upstream of th e 
shock (e.g. lAharonian & Atovan (l_996). lGabici et all (120081) . 
iLee et al1(120081) and lFuiita et alJ (120091) ). However a major¬ 
ity of these particles come from the high-energy portion of the 
7 -ray spectrum and the observation of “low” energy 7 -rays 
may lead to inconsistencies with this scenario. 

5.2.2. Leptonic emission of the observed y-rays 

For inverse Compton scattering to be the dominant mech¬ 
anism producing the 7 -rays of MSH 11-61A, we would re¬ 
quire greater than the entire kinetic explosion energy just in 
electrons, assuming that the electron to proton ratio is similar 
to that measured at Earth (K ep ~ 0.01) and that this emission 
arises from a non-thermal population of electrons being accel¬ 
erated by the shock-front. This makes it difficult to conclude 
that IC scattering is the dominant mechanism producing the 
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Table 4 

Model parameters and density estimates for the pion decay and leptonic model for MSH 11-61A. 


Object 

Distance 

(kpc) 

&proton 

proion 

&elec ^0 

(GeV) 

T?elec 

(GeV) 

Magnetic field 

w 

Ambient density 
«o (cm -3 ) 

X-ray density 
n (cm -3 ) 

MSH 11-61A 

7.00 

4.39 

3.15 6.05 

6.05 

25 

9.20 

see Table [3] 



Figure 8. Broadband fits to the non-thermal radio emissio n (g reen data points), non thermal X-ray upper limit derived from the Suzaku data (dark green limit) 
and the Fermi -LAT 7 -ray emission (as described in Figure[3) of MSH 11-61A. The pion decay, non-thermal bremsstrahlung and synchrotron models defined 
by the parameters in Tableware shown as the solid magenta line, dashed orange line and dot-dashed purple line, respectively. The corresponding IC model falls 
below the plotted axis. 


observed 7 -rays. 

For non-thermal bremsstrahlung to dominate the GeV emis¬ 
sion we require a K ep > 0.2, assuming the maximum density 
derived from our X-ray anal ysis (see Table |3}. Local mea¬ 
surements imply K ep ~ 0.01 dGaisser et al.l 19981) . while 7 -ray 
mode lling of other SNRs imply K ep < 0.01 (e.g. lEllison et al.l 
120101 ) . making it unlikely that non-thermal bremsstrahlung 
emission is the dominant emission mechanism. 

6 . CONCLUSION 

70 months of Fermi-LAT data reveal significant (~ 5tr) 7 - 
ray emission from SNR MSH 11-61A. This emission is con¬ 
sistent with being located along or inside the western limb 
of the remnant given the angular resolution of the Fermi- 
LAT and is adjacent to regions that show a strong recom¬ 
bining plasma component. By modelling the broadband 
spectrum, we find that the emission is best described by a 
hadronic model, while a leptonic scenario is energetically un¬ 
favourable. This is consistent with CO and Hi observations 
that indicate the SNR is interacting with a molecular cloud 
towards the north and southwest. Similar to previous studies, 
the inferred density from our pion decay model is much higher 
than that implied by the thermal X-ray emission. Suzaku data 
reveal that the bulk of the X-ray emission of MSH 11-61A 
arises from a single recombining plasma with enhanced abun¬ 
dances of Mg, Si and S with some regions also requiring an 
underabundance of Ne and Fe, while the emission towards the 
east of the remnant arises from an ionising plasma with Mg, 
Si and S. The origin of the recombining plasma is most like 
adiabatic cooling. We find that the results from our central 
regions (1-3) and our regions 9-10, agree well with those 


that lKamitsukasa et all (120151) obtained for their correspond¬ 
ing regions. The enhancement of Mg, Si and S suggests that 
some of the observed emission arises from shocked ejecta and 
that the progenitor of MSH 11-61A had a mass > 25 M@. 
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